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Abstract 



To realize non-Abelian s-wave topological superfiuid of cold Fermi gases, 
generally a Zeeman magnetic field larger than superfiuid pairing gap is 
necessary. In this paper we find that using an anisotropic and spin- 
dependent optical lattice to trap gases, a new non-Abelian topolog- 
ical superfiuid phase appears, in contrast to an isotropic and spin- 
independent optical lattice. A characteristic of this new non-Abelian 
topological superfiuid is that Zeeman magnetic field can be smaller than 
the superfiuid pairing gap. By self-consistently solving pairing gap equa- 
tion and considering the competition against normal state and phase 
separation, this new phase is also stable. Thus an anisotropic and spin- 
dependent optical lattice supplies a convenient route to realize topo- 
logical superfiuid. We also investigate edge states and the effects of a 
harmonic trap potential. 



PACS number (s): 67.85.Lm, 03.65. Vf, 74.20.-z 

1 Introduction 

In two dimensional condensed matter physics, the search of non-Abelian Majorana fermions 
(MF) [TJ, which are their own antiparticles and may have potential important application 
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in fault-tolerant topological quantum computation [2], attracts much recent attention to 
topological superfluid (TS) and superconductors (TSC) [21 E], which has a full pairing gap 
in the bulk and topologically protected gapless states on the boundary. Inspired by the 
fact that MFs exist at the vortex cores of a two-dimensional (2D) p x + ip y superconductor 
[3], theoretically some practical systems have been proposed, such as fractional quantum 
Hall effect at filling v = 5/2 [5], superfluid He-3 [6], non-centrosymmetric superconductors 
[7J [8], the surface of a three-dimensional topological insulator in proximity to an s-wave 
superconductor [9], [10] and a spin-orbit-coupled (SOC) semiconductor quantum well coupled 
to an s-wave superconductivity and a ferromagnetic insulator [llj . Especially the latter 
has been improved by applying an in-plane magnetic field to a (llO)-grown semiconductor 
coupled only to an s-wave superconductor. This improvement not only enhances the ability 
to tune into the topological superconducting state but also avoid orbital effect of the external 
magnetic field [12]. 

It is widely known that cold Fermi gases can be used to simulate many other systems 
owing to their many controllable advantages and operabilities. Certainly the simulations 
to TS are also possible and as far as what we know is concerned, three main routes have 
been suggested. The first was direct and based on p-wave superfluidity of degenerate Fermi 
gases by means of p-wave Feshbach resonance |13j . Although this method is very simple, it 
is challenging due to short lifetimes of the p-wave pairs and molecules. Subsequently Zhang 
et al. |14j propose to creat an effective p x + ip y TS from an s-wave interaction making 
use of an artificially generated SOC. In fact, SOC have been realized in a neutral atomic 
Bose-Einstein condensate (BEC) and the same technique is also feasible for cold Fermi gases 
[15\ I16j . Realizing that in a dual transformation SOC is formally equivalent to a p-wave 
superfluid gap, Sato et al. [T7] suggest to artificially generate the vortices of SOC by using 
lasers carrying orbital angular momentum. In terms of the latter two ansatzs, in order to 
enter into non-Abelian TS, a Zeeman magnetic field larger than superfluid pairing gap is 
essentially needed. 

In this paper we propose a model to realize non-Abelian TS in cold Fermi gases by 
using an anisotropic and spin-dependent optical lattice (ASDOL) to engineer mismatched 
Fermi surfaces for each hyperfine species. Such two requirements for optical lattices are 
both accessible in an experiment. The anisotropy is determined by the intensity of the 
corresponding pair of laser beams in different directions, while spin-dependence is also 
available in the light of the fact that the strength of the optical potential crucially depends 
on the atomic dipole moment between the internal states involved [18\ [T9l [20l 121] . In 
contrast to an isotropic and spin-independent optical lattice (ISIOL) \14\ 117] . our model 
realizes another new non-Abelian TS phase in which Zeeman magnetic field can be smaller 
than superfluid pairing gap. 

The organization of this paper is as follows. In section 2, we firstly give our model 
and analyze the condition of gap closing. Then we calculate TKNN number Itknn |22] 
of occupied bands addressing the topological properties of the model to obtain topological 
phase diagram. A new non-Abelian TS is discovered. In addition we also investigate the 
properties of edge states to prove bulk-boundary correspondence. In section 3, by taking the 
self-consistency constraint on the s-wave pairing gap into account, the stability of TS and 
the effects of a harmonic trap potential are discussed. We find that this new non-Abelian 
TS is stable. A brief conclusion is given in Section 4. 
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2 Hamiltonian and Topological Phase Diagram 



In this paper we consider s-wave superfluid of cold Fermi gases with Rashba SOC in a two- 
dimension square optical lattice which is anisotropic and spin-dependent. For convenience 
below we assume the lattice constant to be unit. The Hamiltonian suggested is 

H = Y^^ k °- ~ V- ^]% a a k(T + ^2[Jkal t a ki - Aa.^afct + H.C.), (1) 

ka k 

where at creates a spin a =t, i fermion at momentum k = (k x ,k y ). J k = 2J[sm(k y ) + 
is'm(k x )] with J (J > 0) denoting the strength of Rashba SOC. /x, T, A are chemical 
potential, effective Zeeman magnetic field and s-wave superfluid pairing gap, respectively. 

The kinetic energy terms e ka come from anisotropic and spin-dependent optical lattices 
|23| . Imagine tuning the intensities of lasers so that one spin state prefers to hop along the 
x axis and the other prefers to hop along the y axis. In this paper we concentrate on a 
specific situation where the Fermi surfaces of the two spin states are rotated by 90° with 
respect to one another, and only consider a near-neighbor hopping Hamiltonian with single 
particle dispersions 

e fct = -2t a cos (k x ) - 2t b cos (k y ), 

e kir = -2t b cos (k x ) -2t a cos (k y ). (2) 

When there are no Zeeman magnetic field and SOC T = J = 0, this model realized a 
stable paired superfluid state with coexisting pockets of momentum space with gapless un- 
paired fermions [23], similar to the Sarma state in polarized mixtures |24t l25| , [26]. Moreover 
in the strong coupling limit, a d-wave pairing superfluid as well as a d-wave density wave 
state, are also proposed to be achievable in this system [27] . 

Generally speaking, the close of bulk gap is a signal of topological phase transitions, 
although it is not a sufficient condition [28]. Thus to obtain topological phase diagram of 
Hamiltonian (TjQ), we firstly calculate the bulk spectrum of the system to search parameter 
regions for which different topological phases are possible. Then for every regions we cal- 
culate topological invariants to label topological properties. To obtain the bulk spectrum, 
we rewrite the Hamiltonian as 

H = lj24 M ^k, (3) 

1 k 

where ty\ = (aL-, a_^, (i-ki) is a row vector and Mj X 4 is a matrix with M\\ = — M33 = 
efct - 11 - T, M 22 = -M44 = e K - 11 + T, M 12 = M 43 = M 2 *! = M 3 * 4 = J k , M 23 = M 32 = 
-M u = -M41 = A, M13 = M31 = M 24 = M 42 = 0. 

Diagonalizing the matrix M^x^, we find the energy spectrum 

Et = ^/e k+ +e k - + \Jk\ 2 + & 2 ±2E (4) 

with gfc + = -(t a + t b )[cos (k x ) + co s (ky)] - fj,, = (-t a + t b )[cos (k x ) - cos (ky)] - T and 

E, = ^e k+ e k .+e k+ w+e k -^. 
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The close of the bulk energy gap is possible only if E, = 0, in other words 

e k+ +e k - + \Jk\ 2 +& 2 = 2E Q , (5) 

which is equivalent to 

e k+ - fk- - \Jk\ 2 + A 2 = 0, |J fc | 2 A 2 = 0. (6) 

For the s-wave pairing, A / and the second equation in (|6|) is satisfied only when 
k = (0,0), (0,7r), (vr,0), (vr,7r). Substituting these values into the first equation in ([U]), four 
different gap closing conditions are obtained 

ii 2 + a 2 = [2(t b - t a ) + r] 2 , i? + a 2 = [2{t b - t a ) - r] 2 , 

r 2 = A 2 + [2(t a + t b ) + fi] 2 , r 2 = A 2 + [2(t a + t b )-fi} 2 . (7) 

By means of these conditions, Fig. 1(a) shows that there are at least 12 regions, which may 
be topologically distinct. By now we have completed the first step of deciding topological 
phase diagram. Below we will explore the topological numbers to classify the topological 
phases of the model (P). 

In terms of our system, it explicitly breaks the time-reversal symmetry even though 
Zeeman magnetic field is zero. Thus TKNN number Itknn plays a central role in topo- 
logical nature of the system [22j, which is consistent with the conclusion in the periodic 
table of TSC that BdG Hamiltonian only with particle- hole symmetry in two dimension 
is classified by integer group Z [29]. Furthermore the topological properties of the system 
are identified as follows [30]. If TKNN number is nonzero and even (odd), the system is 
Abelian (non-Abelian) TS without (with) non-Abelian anyons; If TKNN number is zero, 
the system is topologically trivial and is normal superfluid. 

TKNN number is defined by 

Itknn = ^~ J d 2 kF l2 (k), (8) 

where Berry connection Ai(k) (i = 1,2) and associated field strength F\ 2 {k) are given by 

Mk) = E < <f>n(k)\d ki <j>n(k) > 
E n <0 

F 12 (k) = d l A 2 (k) - d 2 A x {k) (9) 

with \4> n (k) > being a normalized wave function of the nth band such that M^x4,\(f>n{k) >= 
E n (k)\4> n (k) >. It is should be noted that the sum in ([9]) is restricted to occupied bands. 

TKNN numbers of 12 regions have been numerically calculated [31J and given in Fig. 1(a). 
Totally 5 non-Abelian TS and one Abelian TS phases are found. In order to make a contrast, 
in Fig. 1(b) we also show the topological phase diagram of an ISIOL and there are 4 non- 
Abelian TS and 1 Abelian TS phases. By comparison the effects of ASDOL are mainly 
twofolds. On the one hand it creates another non-Abelian TS phase around the chemical 
potential /i = 0, which can be realized for Zeeman magnetic field T less than superfluid 
pairing gap A. This is our main result in this paper and can be understood from the 
fact that an ASDOL (t a ^ t b ) supplies an effective Zeeman magnetic field ±2(t a — t b ), as 
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seen from (7). On the other hand in the direction of increasing Zeeman magnetic field it 
separates two successive non-Abelian TS phases in Fig. 1(b). 

According to the bulk-boundary correspondence, a topologically nontrivial bulk guar- 
antees the existence of topologically stable gapless edge states on the boundary. Cold Fermi 
gases with sharp edges may be realized along the lines proposed in [32]. In our case, the 
gapless edge states is a chiral Majorana fermion mode. It should also be remembered that 
the core of a vortex is topologically equivalent to an edge which has been closed on itself. 
The topologically protected edge modes we describe here are therefore equivalent to the 
Majorana fermions known to exist in the core of vortices of p-wave superfluid [33] . 

To study edge states, we transform Hamiltonian ([1]) into lattice representation 

H = -t a ^2 [ a !t°*+zt + a li a i+yl + H.C.] -tb^2 [a^cii+yf + al^a i+x ± + H.C.] 

i i 

+ J X^( a It a *+z| ~ 4l a i+xt) ~ K a ^ a i+yi + a ll a i+yt) + H.C.] 

i 

- + aT)4 a a ia - Aj2(4 t 4i + B.C.) (10) 

ia i 

where a\ a is the creation operator of a fermion with spin a at lattice site i = (i x ,iy)- 
Without loss of generality we suppose that the system has two open boundary in the re- 
direction and is periodic in the y-direction. After performing a Fourier transformation along 
the y-direction only, we numerically diagonalize the Hamiltonian (|1U|) and correspondingly 
obtain the excitation spectrum E n ^ y with subscript n labeling different energy levels. 

In Fig. 2 we shows the energy spectra for an ASIOL with edges at i x = and i x = 30. 
These 12 figures correspond to 12 regions in Fig.l(a), respectively. Fig.2(b), (d), (f), (h), 
(j) are for non-Abelian TS phases, (k) for Abelian TS phase and others non-TS phases. 
As stated in [30, 33J, we also find that in the non-Abelian TS phases, where Itknn = ±1 
and (— lyrKNN — _i ) a s i n gl e p a i r G f gapless edge modes appears. Since on a given 
boundary there is no state available for backward spinconserving scattering, these states 
are topologically protected [33]. In this case the edge zero mode is a chiral Majorana 
fermion. While in Abelian TS phase Itknn = 2 and non-TS phases Itknn = 0, where 
(— lyrKNN — ]_ ) the gygtem contains either zero or two pairs of edge states. These edge states 
are topologically trivial for either edge modes can perform spin-conserving backscattering. 
In addition it is also found that edge excitations only cross at zero energy with a linear 
dispersion at k y = or tt. Their origin is topological and can be identified by a topological 
winding number I{k y ) defined only for k y = and tt [30]: when I{k y ) is non-zero for k y = 
or 7T, the energy of the gapless edge mode becomes zero at this value of k y . 

3 Phase Diagram with Self-Consistent Pairing Gap 

In section 2 we have determined the topological phase diagram for a fixed s-wave superfluid 
pairing gap. However for a realistic physical system the pairing gap generally varies when 
the chemical potential and Zeeman magnetic field change. Especially Zeeman magnetic field 
breaks time-reversal symmetry and weakens the stability of superfluidity. A well-known ex- 
ample is so-called Chandrasekar-Clogston (CC) limit [3H [35] in superconducting systems 
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without SOC. Hence it is possible that not all phases in Fig.l(a) are accessible. In this sec- 
tion within BCS mean-field theory, we self-consistently determine s-wave superfluid pairing 
gap and consider the competition from normal phase and phase separation to investigate 
the stability of TS. 

Let — U (U > 0) denote the effective attraction strength between fermions, then the pair- 
ing gap A = UJ2k < a -kl a k j [ > can be obtained from the minimization of thermodynamic 
potential Q s = J2k £fc+ — \i-^t + ^k) + NA 2 /U . The unstability of superfluidity against 
phase separation is signalled by the condition A / and d 2 Q s /dA 2 < 0, while unstability 
against normal state is Q n < Q s , where U n is thermodynamic potential of normal state. For 
parameter chosen superfluidity is robust as long as A / 0. The resulting zero-temperature 
pairing gap and phase diagram are shown in Fig. 3. Fig. 3(b) shows similar structures with 
Fig. 1(a) except TS with large Zeeman magnetic field are not available. Comparing Fig.3(a) 
with (b), it is easily seen that for small chemical potential and Zeeman magnetic field, 
although the pairing gap has a larger magnitude (A ~ t a ) than Zeeman magnetic field, 
the system can be non-Abelian TS with Itknn = — 1- This finding is very important and 
suggests that such a new non-Abelian TS we have found in section 2 is stable. 

At last we consider the effects of a harmonic trap potential by local density approxi- 
mation (LDA). Under this approximation the system is locally uniform and trap potential 
effectively provides a scan over the chemical potential, in other words local chemical poten- 
tial [1r = fi — AtR 2 with At denoting the strength of an isotropic potential. From Fig. 3(b) 
depending on Zeeman magnetic field some different topological phases can coexist in dif- 
ferent concentric spherical shells. In Fig.4 we show the space distributions of pairing gap, 
particle number and TKNN number. It is should be noted that on the boundary of different 
topological phases there always exists zero-energy states due to the change of topological 
invariant. 



4 Conclusions 

In conclusion we have investigated the effects of an anisotropic and spin-dependent optical 
lattice (ASDOL) on non-Abelian topological superfluid and found that a new non-Abelian 
topological superfluid phase exists steadily, in contrast to an isotropic and spin-independent 
optical lattice. Moreover in this new non-Abelian topological superfluid phase Zeeman 
magnetic field can be smaller than the superfluid pairing gap. Thus an ASDOL supplies a 
convenient route to realize TS. In addition we also calculated chiral Majorana edge states 
and investigated the effects of a harmonic trap potential. 
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Figure 1: Topological phase diagram. The numbers in different regions are TKNN numbers. 
We have chosen tf,/t a = 0.35, J/t a = 0.5 and A/t a = 0.8. (a) is for an anisotropic and spin- 
dependent optical lattice and (b) is for an isotropic and spin-independent optical lattice. 
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Figure 2: The energy spectra of an anisotropic and spin-dependent optical lattice with open 
edges at i x = and i x = 30. We have chosen tb/t a = 0.35, J/t a = 0.5 and A/t a = 0.8. 
(a) n/ta = -4.0, T/t a = 1.0, (b) n/t a = -4.0, T/t a = 2.5, (c) n/t a = -4.0, T/t a = 4.0, 
(d) M /t a = -4.0, r/t a = 6.5, (e) p/t a = 0.0, r/t a = 6.5, (f) /u/t a = 4.0, T/t a = 6.0, (g) 
li/t a = 4.0, r/t a = 4.0, (h) n/t a = 4.0, T/t a = 2.0, (i) /u/t a = 4.0, T/t a = 1.0, (j) fjt/t a = 0.0, 
r/t« = 1.0, (k) fj,/t a = 0.0, T/ta = 2.5, (1) fi/t a = 0.0, r/i a = 0.2. These figures correspond 
to different regions in Fig. 1(a), respectively. 
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Figure 3: s-wave superfluid pairing gap (a) and topological phase diagram (b) from self- 
consistent mean-field solution. In (b) grey, black, blue and purple colors correspond to 
Itknn = 0, —1, 1, —2 respectively. We have chosen tb/t a = 0.35, J/t a = 0.5 and U/t a = 4.0. 
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Figure 4: The space distributions of pairing gap, particle number and TKNN number. We 
have chosen t b /t a = 0.35, J/t a = 0.5, [i/t a = 3, A t /t a = 0.01, T/t a = 1.2 and U/t a = 4.0. 
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